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1. The strange sigma term 

To quantify and investigate the nature of chiral symmetry breaking in QCD one would ideally 
like to measure objects such as the matrix element of the symmetry-breaking Hamiltonian density 
between nucleon states. While this is not possible directly, the sigma term can be obtained from 
on-shell scattering amplitudes. 

The baryonic sigma terms are defined as scalar form factors: for each quark flavor q and 
baryon B, 



For the nucleon, the so-called kN sigma commutator and the strange sigma commutator are defined 

as 



where mi = (m u + nid)/2. 

Written in this form, it is clear that the sigma terms measure the strength of the various matrix 
elements rriqq {q = u,d,s) in the baryons. Hence, these terms quantify both the contribution of 
chiral symmetry breaking to the baryon masses, and the scalar density of quarks in the baryons. 
As many dark matter candidates have interactions with hadronic matter which are determined by 
couplings to the light and strange quark sigma terms, these commutators are also of essential im- 
portance to the interpretation of experimental searches for a particle candidate of dark matter [1-6]. 

Largely because of this importance to dark matter searches, recent years have seen a concerted 
effort towards precisely determining the nucleon sigma terms in particular. While the best value for 
o k n has remained somewhat stable - modern lattice calculations support a value « 45 MeV 
as determined from nN scattering through a dispersion relation analysis [7,8] - the best value for 
a s has seen an enormous revision. Traditional indirect evaluations using o„n and a best-estimate 
for the non-singlet contribution Go = mi(N\uu + dd — 2ss\N) have yielded a value for o s as large 
as 300 MeV. Since the strange sigma commutator may be interpreted as the contribution to the 
mass of the nucleon from the strange quark, a value as large as this would indeed be remarkable; 
it would suggest that almost a third of the nucleon mass arises from non-valence quarks. This 
appears incompatible with the widely used constituent quark models, for example. While early 
lattice simulations appeared to support large values of o s , advances in lattice QCD have revealed a 
strange sigma term of 20-50 MeV [9, 1 1-16,26], an order of magnitude smaller than was previously 
believed, and significantly more precise. Recent work from QCDSF supports the suggestion that 
operator mixing is responsible for erroneously large values of o s in the early lattice studies [11]. 

Shanahan et al. [17] recently presented new values for the octet baryon sigma terms. These 
were extracted using the Feynman-Hellman theorem from a chiral perturbation theory fit to octet 
baryon mass data from the PACS-CS collaboration [22]. Of particular interest is the precise deter- 
mination of the strangeness sigma term o s = 21 ±6 MeV. In this work, we elaborate on the critical 
effect which the choice of scale setting has on this value. 



C Bq = m q {B\qq\B). 



(1.1) 



Gjin = mi{N\uu + dd \N) 
G s =m s (N\ss\N), 



(1.2) 
(1.3) 
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1.1 Scale setting and the Feynman-Hellman theorem 

A common feature of lattice simulations is that all quantities are calculated in units of the 
unknown lattice spacing a, which must be determined by matching an observable to its physical 
value. This can be done in a variety of ways. Two common methods, often referred to as 'mass 
independent' and 'mass dependent' scale setting schemes, are of interest to us here: 

1. The inverse bare coupling j8 determines the lattice spacing a. That is, points at fixed j3 
are extrapolated to the physical point (usually linearly in the bare quark mass am qse J, and 
the value of some observable at that point is used to set the common scale a for all lattice 
ensembles at that common /3 . 

2. The lattice spacing varies with bare quark mass. That is, the lattice spacing a is determined 
separately for each set of bare parameters (j8 , am qscii ) by using a physical observable that is 
assumed to be independent of the quark masses (often the Sommer scale ro). 

We could think of these two choices as different ways of absorbing the observed quark mass 
dependence of the ratio r$/a at fixed /3. Method 1 essentially assumes that this dependence may 
be attributed to the variation of ro with quark mass, while method 2 instead assumes that a is 
changing. What is not often considered is that both ro and a may have some dependence on the sea 
quark mass, which would lead to an intermediate scale in some sense. Such a 'mixed' scale setting 
procedure would be non-trivial to implement. 

Of course, in the continuum limit, and when the chiral extrapolation to physical quark masses 
has been performed, the results of each method of scale setting must agree for physical observables. 
When considering derived quantities which are expressed as derivatives with respect to quark mass, 
however, the choice of scale setting method becomes far more significant; these quantities, by the 
very definition of the derivative, depend on the scale away from the physical point and hence on 
the scale setting scheme chosen. The Feynman-Hellman theorem expresses the baryon sigma terms 
precisely as such derivatives [23], 

dM B 



o Bq = m q - — . (1.4) 

om q 



Clearly there are two primary methods of evaluating the baryon sigma terms on the lattice. The 
first is by direct calculation of the relevant matrix elements. This method does not suffer from any 
major scale setting ambiguities. The other is to use lattice data for the octet baryon mass spectrum 
and the Feynman-Hellman theorem to determine the sigma terms - the 'spectrum' method. The 
results of this method will depend on the choice of scale setting procedure, as discussed above. 

This effect may potentially be significant for the strange sigma term in particular. One reason 
is the well-known 'beta shift' effect which is observed when unquenching lattice simulations; at 
fixed /3, the ratio ro/a increases when dynamical quarks are added [24], which can be interpreted as 
a sea quark dependence of either ro, or the spacing a. As this shift can be significant, the choice of 
scale setting absorbs a possibly large effect, and hence will lead to non-negligible differences in the 
results of derivatives with respect to sea quark mass calculated with each of the two choices. Addi- 



tionally, the strange quark is considerably heavier than the light quarks, which, by Equation (1.4), 
serves to amplify the effect which a change in the slope of baryon mass with respect to strange 
quark mass has on the sigma term. 
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To illustrate the significance of the choice of scale setting method, we consider the recent 
determination of the octet baryon sigma terms presented by Shanahan et al. [17]. This work used 
method 2 to set the scale. If an identical analysis is repeated, using the same PACS-CS lattice 
data but with the scale set using method 1 instead, the results are significantly different. While the 
chiral extrapolations of the octet baryon masses are compatible at the physical point, as expected, 
and the qualities of the two fits are similar (# 2 /d.o.f 0.41 and 0.67 for methods 2 and 1), the value 
of the strangeness sigma term changes significantly. In particular, we note that the value for a s 
in the nucleon changes from 21 ±6 MeV to 59±7 MeV as the scale setting scheme is changed 
from method 2 to method 1. Despite the small statistical and model-dependent uncertainties of 
the calculation, there is clearly a significant dependence on the scale setting method. Results for 
the light quark sigma terms and for other octet baryons are given in tables [I] and ||. For each 
baryon, there is a shift of about 35 MeV in G Bs . Clearly, a deeper understanding of scale setting on 
the lattice is of essential importance for the interpretation of spectral determinations of the sigma 
terms. 



B Mass (GeV) Experimental a B i (GeV) o Bs (GeV) 

N 0.959(24)(9) 0.939 0.045(5)(5) 0.021(6)(0) 

A 1.129(15)(6) 1.116 0.029(3)(2) 0.159(7)(1) 

I 1.188(H)(6) 1.193 0.024(2)(2) 0.204(8)(1) 

S 1.325(6X2) 1-318 0.0117(9)(5) 0.316(9)(2) 

Table 1: Extracted masses and sigma terms for the physical baryons, with the lattice scale set using method 
2. The first uncertainty quoted is statistical, while the second results from the variation of various chiral 
parameters and the form of the UV regulator as described in [17], from which these results are taken. The 
experimental masses are shown for comparison. 



B Mass (GeV) Experimental a B i (GeV) o Bs (GeV) 

N 0.944(22)(9) 0.939 0.051(3)(6) 0.059(6)(1) 

A 1.119(15)(6) 1.116 0.034(2)(2) 0.193(8)(1) 

I 1.182(H)(6) 1.193 0.028(2)(2) 0.241(9)(2) 

S 1.323(7X2) 1.318 0.0151(9)(4) 0.352(H)(1) 

Table 2: Extracted masses and sigma terms for the physical baryons, with the scale set using method 1. The 
uncertainties are as in table [j]. 



The effect of the choice of scale setting on such spectral determinations of the sigma terms has 
been discussed in the literature. Both [18] and [19] argue for method 1, in the former case based 
on observations of scaling violation, while others, for example [20, 21] instead favour method 2. 
There is no clear consensus as to the most appropriate way to set the scale. 

In the work referenced here [17], Shanahan et al. chose method 2 to set the scale for their 
primary analysis, using the argument that the Sommer scale tq, based on the force between static 
quarks at some distance, should be independent of small changes in the (small) bare quark masses. 
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This approach is strongly supported by the fact that the results of this calculation agree with direct 
lattice calculations at unphysical quark masses. For example, as discussed in [17], the QCDSF 
Collaboration [11] recently presented a value ct Sqcdsf = 12+jg MeV at (m n ,m.K) = (281,547) MeV 
which agrees well with the results of the calculation of [17] with the scale set using method 2, 
namely a« 2 = 16(5)(1) MeV. In contrast, the direct result is not consistent with the value o Sl = 
56(6) (1) MeV obtained when the scale is set using method 1. Similarly, other recent direct lattice 
calculations favour smaller values of a s [25, 26]. While [15] present a larger o s , naively more 
compatible with scale setting method 1, the authors point out that the effects of operator mixing 
likely serves to increase their result, again suggesting a true value more compatible with method 2. 
We point out that while most recent direct results present small values for a s , not all calculations 
do support this trend [27,28]. 

Despite this debate, we suggest that method 2 of setting the lattice scale yields results which 
are more consistent with direct calculations. Of course, we cannot rule out the possibility of 'mixed 
scale setting' as discussed earlier, and the precision of the result of a s = 21 ± 6 MeV relies heavily 
on the assumption that ro does not change with quark mass. We emphasize that as more precise 
direct lattice calculations for the strangeness sigma commutator appear, we will not only more 
precisely determine the value of this term as needed for dark matter calculations, but we stand to 
gain significant insight into the problem of scale setting on the lattice and indeed into QCD itself. 

2. Up, Down and Strange Sigma Terms 

As discussed in previous sections, Shanahan et al. [17] recently presented new values for the 
light and strange quark octet baryon sigma terms. This work applied the Feynman-Hellmann re- 
lation to the results of a chiral extrapolation of PACS-CS Collaboration [22] lattice data for octet 
baryon masses. Arguably the most significant result presented was a small value for the strangeness 
sigma term, namely o s = 21 ± 6 MeV at the physical point. In what follows, we summarize addi- 
tional results of this work which were not presented in the original paper. In particular, using the 
approach presented in [29], individual up, down and strange sigma terms as relevant to supersym- 
metric dark matter searches [1, 33, 34] may be evaluated. 

In reference [29], the authors used the same chiral perturbation theory fits to octet baryon mass 
data as was used in their evaluation of the sigma terms [17]. By extending the chiral extrapolation 
to the case of broken isospin symmetry, i.e., m u ^ rtid, they determined the mass splittings between 
members of baryonic isospin multiplets based on their fit to isospin-averaged lattice data. We use 
precisely the same method to extract the individual, isospin-broken, sigma terms fj u / fjd- Just as 
was done in [29], we take as input the physical up-down quark mass ratio [30], 

R:= — =0.553 ±0.043, (2.1) 

determined by a fit to meson decay rates. We note that this value is compatible with more recent 
estimates of the ratio from 2 + 1 and 3 flavor QCD and QED [31, 32]. All calculations include in 
quadrature the uncertainty due to the stated range for R. 

Following the notation of Ellis et al. [1, 33, 34], we define individual quark (q) sigma terms in 
the proton as 

m pfTq = {p\ m qW\P) = m q B qi ( 2 - 2 ) 
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(2.3) 



These terms are related to the previously defined o s and o„n by 



a KN = l -(m u + m d )(B u P) +*S°) 



(2.4) 



(2.5) 



The results of our analysis are 




0.019(3), 
0.015(2), 



0.027(4), 
/g= 0.033(5), 




An) 



0.023(7), 
0.022(6). 



(2.6) 



(2.7) 



We also find 



z= 1.27(3) 



« 0.79(2), 



jj-{N\uu - dd\N) = 0.0049(2). 



(2.8) 



The quoted errors include correlated uncertainties between all fit parameters, and allow for some 
variation of phenomenologically set quantities including the up-down quark mass ratio R. For 
details see [17,29]. 

These results may be compared with those in [34] (we note that the same R value was used in 
that analysis). While our up and down sigma terms are compatible, within uncertainties, with these 
results, we point out that the strange sigma terms fj s and the parameter z resulting from our work 
are significantly smaller than the previous values. 

Of course, the discussion of the previous section as regards a s applies equally to f Ts : setting 
the lattice scale with fixed coupling j3 corresponding to a fixed lattice spacing yields far larger 
values for the strange sigma terms, namely f T ^ = 0.063(7) and f T n J = 0.062(7), but a compatible 
value z = 1.25(2). 
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